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Abstract  —  We  conjecture  that  good  column’hased  place¬ 
ments  can  be  produced  by  minimizing  two  wire  crossing  num¬ 
bers:  (1)  the  total  wire  crossing  of  all  edges  between  cells  in 
the  wiring  channel j  and  (2)  the  maximum  wire  crossing  on  the 
imaginary  outlines  that  separate  cells  on  the  opposite  ends  of 
the  channels.  We  leverage  the  canonical  form  of  the  multi¬ 
level  bipartite  directed  graph  to  formalize  a  unit-grid  model 
that  allows  us  to  define  and  evalutate  parameters  such  as  total 
wire  crossing,  critical  wire  crossing,  total  wire  length,  critical 
wire  length,  critical  wire  density,  total  wire  density,  as  well  as 
height,  width,  and  area  of  the  embedded  graph. 

We  implemented  a  prototype  placement  algorithm  TOCO  that 
minimizes  the  cost  of  wire  crossing,  and  a  universal  unit-grid 
based  placement  evaluator  place^eval.  We  have  designed  a 
number  of  statistical  experiments  to  demonstrate  the  feasibility 
and  the  promise  of  the  proposed  approach. 

Keywords;  wire  crossing,  wire  length,  cell  placement,  bipar¬ 
tite  graph  model,  design  of  experiments,  benchmarking. 

L  Introduction 

Area  and  performance  in  submicron-technology  VLSI  circuits 
are  critically  dominated  by  the  interconnect.  The  variations 
in  interconnect  that  arise  after  cell  placement  and  routing, 
such  as  wire  crossings,  wire  density,  wire  loading,  wire  length, 
etc.,  may  be  hard  if  not  impossible  to  predict.  Similarly,  the 
interconnect  models  associated  with  the  logic  implementa¬ 
tion  and  logic  testability  prior  to  placement  and  routing  may 
change  considerably  after  cell  placement  and  routing. 

In  row-based  cell  placement,  the  cost  function  that  con¬ 
tinues  to  dominate  the  research  relates  mostly  to  total  wire- 
length,  e.g.  the  simulated  annealing  approach  [1].  The  length 
of  the  interconnect  in  randomized  and  optimized  placements 
has  been  formalized  in  [2].  Abstract  placement  models,  such 
as  presented  in  [3],  also  aim  to  minimize  the  channel  density. 

The  crossing  number  and  wire  area  have  been  investigated 
for  effective  wire  lower  bounds  and  edge  length  in  a  variety  of 
computational  VLSI  circuits,  using  a  graph-based  grid  model 

[4] .  In  this  paper,  we  propose  a  new  unit-grid  model,  de¬ 
rived  from  the  a  multi-level  bipartite  graph  in  canonical  form 

[5] .  We  conjecture  that  good  placements  can  be  produced  by 
minimizing  two  wire  crossing  numbers  in  our  bipartite  graph 
model:  (1)  the  total  wire  crossing  of  all  edges,  and  (2)  the 
maximum  wire  crossing  of  all  cell  channel  sepctrator  segments. 

Crossing  theory  has  been  developed  to  improve  the  read¬ 
ability  of  hierarchical  structures  [6].  The  problem  for  placing 
the  nodes  for  minimum  wire  crossing  is  NP-complete  [7],  even 
for  2-layer  graphs.  Several  heuristics  -  barycenter  [8],  median 
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[9],  assignment  [10]  -  have  been  proposed  to  minimize  the 
wire  crossings.  A  survey  of  exact  and  heuristic  algorithms  is 
available  in  [11].  An  improvement  on  the  median  heuristic  has 
been  described  in  [12]  and  we  use  the  associated  drawing  pro¬ 
gram  DDT  to  support  our  initial  experiments  reported  in  this 
submission.  We  are  also  investigating  alternative  approaches 
to  minimize  wire  crossing,  subject  to  constraints  specific  to 
VLSI  circuits  and  expect  to  report  on  them  in  the  final  sub¬ 
mission. 

We  argue  that  for  layout  problems  in  particular,  case-by- 
case  evaluations,  however  detailed,  of  a  few  unrelated  bench¬ 
mark  circuits  are  not  likely  to  reveal  a  set  of  statistically  con¬ 
sistent  results  that  can  drive  and  support  important  improve¬ 
ments  for  the  new  generation  of  algorithms  and  tools.  To  il¬ 
lustrate  our  point,  we  analyze  grid-based  layouts  as  well  as 
placed- and-routed  layouts  for  a  number  of  equivalence  mutant 
classes  introduced  in  [5]^.  For  convenience  of  the  reviewers, 
we  copy  the  section  on  the  design  of  experiments  and  the  mu¬ 
tant  classification  we  used  in  [5].  In  the  the  final  submission, 
this  section  will  be  edited  to  minimum  size. 

The  paper  is  organized  into  the  following  sections  : 

(2)  background  and  motivation; 

(3)  correlations  in  2-layer  graphs; 

(4)  notation  and  definitions; 

(5)  cell  and  net  folding  and  compaction; 

(6)  design  of  experiments; 

(7)  summary  of  experiments; 

(8)  conclusions. 

IL  Background  and  Motivation 

We  use  a  simple  11-node  planar  circuit  to  informally  intro¬ 
duce  the  proposed  unit-grid  model  and  illustrate  the  conjec¬ 
ture  that  good  placements  can  be  produced  by  minimizing  two 
wire  crossing  numbers  in  its  directed  bipartite  graph  model 
representation.  A  summary  of  our  experiments  is  shown  in 
Figure  1: 

(a)  The  circuit  is  not  drawn  in  the  planar  form.  Rather, 
we  draw  it  as  an  instance  of  a  random  placement,  in  a 
canonical  form,  of  a  bipartite  directed  acyclic  graph  [5]. 
The  ‘square’  and  ‘triangle’  nodes  represent  net  nodes, 
the  circles  and  ellipses  are  cell  nodes.  Each  circle  rep¬ 
resents  a  logic  cell  in  the  original  circuit,  each  ellipse 
represents  a  feedthrough  cell  introduced  by  the  canonical 
form.  We  call  the  horizontal  lines,  drawn  as  dotted  lines 
and  indexed  from  0  to  4,  cell  separators.  We  call  the 
vertical  lines  that  could  be  drawn  through  the  net  nodes 
and  indexed  from  0  to  3,  net  separators.  The  total  wire 
crossing  number  for  all  edges  in  the  graph  as  drawn  here, 
is  54.  Starting  with  the  placement  shown,  not  even  the 
state-of-the-art  algorithm  such  as  DOT  [12]  can  render  the 
circuit  planar  (wire  crossing  is  minimized  to  a  value  of 

7)* 

(b)  Analysis  of  the  left-most  net-to-cell  channel  in  (a) 
shows  a  total  wire  crossing  of  23,  wire  length  of  23.5, 
and  the  critical  wire  density  of  7.  Wire  crossing  refers  to 

^Available  as  a  preprint  from  http://www.cbl.ncsu.edu/publications/). 
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Fig.  1.  Subproblems  in  circuit  layout:  experiments  with  mutant  Class  ‘D’  of  circuit  ‘planar ^ 


crossing  of  edges  between  net  nodes  and  cell  nodes.  If  an 
edge  between  two  adjacent  nodes  crosses  no  cell  separa¬ 
tors  (shown  as  dotted  lines),  we  define  its  length  as  0.5, 
and  increase  its  length  by  1  whenever  it  crosses  a  cell 
separator.  The  critical  wire  density  is  evaluated  as  the 
maximum  number  of  edges  crossing  the  cell  separator, 
(c)  Wire  crossing  minimization  of  net-to-cell  channels  in 
(b)  reduces  total  wire  crossing  to  0,  wire  length  to  8.5, 


and  the  critical  wire  density  to  1. 

(d)  An  instance  of  the  standard  cell  layout  with  OASIS 
[13]  shows  a  number  of  wire  crossings.  In  fact,  upon 
translating  the  specific  OASIS  placement  into  the  bipar¬ 
tite  graph  form  on  the  same  grid,  we  have  1  wire  crossing 
in  channels  (0-1),  3  wire  crossings  in  channels  (1-2),  and 
1  wire  crossing  in  channels  (2-3).  Inspecting  the  actual 
layout,  we  find  more  wire  crossings  in  each  of  the  chan- 
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nels. 

(e)  Table  summarizing  an  experiment  with  100  instances 
of  the  netlist  of  the  11-node  planar  circuit,  each  in  ran¬ 
dom  order,  before  submission  to  wire  crossing  minimiza¬ 
tion  with  DOT  [12].  This  is  the  W-ClassD  circuit  class, 
a  special  class  of  isomorphic  mutants  in  the  context  of 
the  design  of  experiments  as  introduced  in  [5].  We  mea¬ 
sure  the  initial  and  final  value  of  sample  mean,  sample 
variance,  and  95%  confidence  interval  of  the  mean  for 
several  parameters  we  have  already  defined  or  will  de¬ 
fine  later  in  the  paper:  total  wire  length  (TWL),  critical 
wire  length  (CWL),  total  wire  crossing  (TWC),  width  of 
placement  (Width),  height  of  placement  (Height),  and 
area  of  placement  (Area). 

(f )  Histogram  of  circuit  area  for  100  instances  of  the  netlist 
of  the  11-node  planar  circuit,  each  in  random  order,  sub¬ 
mitted  for  layout  optimization  with  OASIS  [13]. 

(g)  Histogram  of  circuit  minimized  wire  crossing  for  100 
instances  of  the  netlist  of  the  11-node  planar  circuit,  each 
in  random  order,  submitted  for  wire  crossing  minimiza¬ 
tion  with  DOT  [12]. 

Conclusions  that  we  draw  from  the  observations  recorded  in 
Figure  1  include: 

(a)  The  circuit  is  embedded  onto  a  number  of  unit-size 
grid  regions.  A  formal  description  is  given  in  Section  IV. 
The  bipartite  subgraph  consisting  of  net  nodes  on  the 
left  and  cell  nodes  on  the  right  is  called  the  net-to-cell 
channel.  Similarly,  the  bipartite  subgraph  consisting  of 
cell  nodes  on  the  left  and  net  nodes  on  the  right  is  called 
the  cell-to-net  channel 

(b)  Random  placement  of  nodes,  either  in  a  net-to-cell  or 
a  cell-to-net  channel,  can  result  in  a  large  number  of  wire 
crossings. 

(c)  Wire  crossing  minimization,  either  in  a  net-to-cell  or  a 
cell-to-net  channel,  can  reduce  not  only  the  wire  length 
but  also  the  wire  density. 

(d)  Traditional  layout  optimization  algorithms  are  not 
minimizing  the  cost  of  wire  crossing. 

(e)  Wire  crossing  minimization  of  a  complete  circuit,  not 
only  a  single  channel  segment,  can  improve  a  number  of 
layout  parameters  in  the  bipartite  graph  model,  not  only 
wire  crossing. 

(f)  Minimized  circuit  area  for  100  instances  of  the  netlist 
of  the  11-node  planar  circuit,  each  listed  in  different  ran¬ 
dom  order,  is  a  random  variable  and  may  be  far  from  best 
possible. 

(g)  Minimized  wire  crossing  for  100  instances  of  the  netlist 
of  the  11-node  planar  circuit,  each  listed  in  different  ran¬ 
dom  order,  is  a  random  variable.  Only  relatively  few 
netlists  have  been  returned  with  the  known  minimum 
wire  crossing  of  0.  Unlike  in  the  case  of  layout  in  (f), 
this  distribution  is  not  normal  since  more  solutions  are 
clustered  to  the  left  of  the  minimum  value  ~  lending  sup¬ 
port  to  our  conjecture  that  optimizing  the  placement  by 
minimizing  wire  crossing  with  a  good  algorithm  may  be 
as  effective,  if  not  more,  than  minimizing  the  wire  length 
estimates. 

Before  we  proceed  with  formal  notation  and  definition,  we 
digress  with  one  more  experiment,  this  time  scaling  the  prob¬ 
lem  for  a  2-layer  graph  minimization  algorithm. 

III.  Correlations  in  2-Layer  Graphs 

The  2-layer  subgraphs  in  the  example  shown  in  Figure  1  are 
too  small  to  analyze  correlations  between  grid-based  variables 
such  as  wire  crossing,  wire  length,  and  wire  density. 

The  net-to-channel  example  in  Figure  1(c)  is  a  special  case 
of  a  two-layer  graphs  from  several  families  of  parameterized 


2-layer  directed  sparse  graphs:  Vn+lE3-n  (E3-graphs)  [14]. 
Properties  of  E3-graphs  are:  number  of  nodes  at  level  0  = 
n+1,  number  of  1-input  nodes  at  level  1  =  3,  number  of  2-input 
nodes  at  level  1  =  n,  number  of  edges  =  2*n  +  3,  number  of 
wire  crossings  =  0. 

We  have  chosen  to  analyze  a  graph  with  65  net  nodes,  131 
edges,  and  3+64  cell  nodes  (3  single-input  nodes,  64  2-input 
nodes).  We  label  its  reference  graph  as  V65E3-64.  Next,  we 
created  8  wiring  signature-invariant  equivalence  classes,  100 
circuits  in  each  class  (each  with  65  net  nodes,  131  edges,  and 
3+64  cell  nodes)  [5]. 

Here,  we  again  analyze  placements  of  the  graph-isomorphic 
class  WXlassD.  Evaluations  of  all  100  such  random  placements 
are  shown  in  Figure  2(a-c):  note  the  near-normal  distribution 
of  wire  crossings  (ranging  from  3400  to  4800)  and  the  near 
perfect  correlation  with  wire  length.  As  to  the  maximum 
wire  density,  it  can  vary  from  59  to  89  wires. 

Next,  we  submit  the  same  100  circuits  to  the  drawing  pro¬ 
gram  DOT  [12].  While  much  improved,  only  14  out  of  100  cir¬ 
cuits  could  be  placed  with  wire  crossing  of  0,  optimum  wire 
length  of  194.5,  and  the  maximum  wire  density  of  1.  The 
remaining  circuits,  while  much  better  placed  than  the  initial 
placement,  are  subopt imal  and  can  still  vary  over  a  relatively 
large  range. 

In  conclusion,  the  example  in  this  section  re-enforces  our 
conjectures  based  on  the  earlier  experiment  in  Figure  1: 

•  It  is  a  fallacy  to  rely  on  a  single  measurement  of  any 
benchmark  circuit  -  variations  for  many  of  the  ‘improve¬ 
ments’  published  to  date  may  well  be  attributed  to 
chance  rather  than  any  intrinsic  improvement  of  the  al¬ 
gorithm. 

•  The  correlation  coefficient  between  the  total  wire  cross¬ 
ing  and  total  wire  length  is  high  before  and  after  wire 
crossing  minimization;  minimizing  total  wire  crossing 
minimizes  total  wire  length. 

•  The  correlation  coefficient  between  the  total  wire  cross¬ 
ing  and  maximum  wire  density  is  much  higher  for  ran¬ 
dom  placement  than  it  is  for  the  near-optimal  placement. 
It  should  be  possible  to  ‘improve’  maximum  wire  den¬ 
sity  without  a  major  change  in  wire  crossing  (and  wire 
length). 

•  Any  improvements  in  (1)  the  total  wire  crossing  of  all 
edges,  and  (2)  the  maximum  wire  crossing  of  all  cell 
channel  separator  segments  may  translate  to  improved 
layout  area  and  performance. 

IV.  Notation  and  Definitions 

The  traditional  graph-based  models  of  a  directed  netlist  are 
not  effective  for  the  problem  we  want  to  consider  in  this  pa¬ 
per.  On  the  other  hand,  a  model  of  a  netlist  as  a  directed 
hypergraph  is  not  unique.  The  important  items  in  this  work 
are  the  notion  of  cell  level,  levels  of  net  pins,  and  netspan  . 
The  canonical  form  of  a  bipartite  directed  graphs  a  multi-level 
graph  structure  of  alternating  sets  of  net  nodes  and  cell  nodes, 
is  a  simple  transformation  of  the  underlying  netlist:  levels  of 
some  of  its  pins  are  redefined,  and  a  new  type  of  cell  node,  a 
feedthrough  cell  is  introduced  [5]. 

Figure  3  relates  to  the  definitions  of  a  grid  region  R{iyj). 
The  vertical  center-line  of  each  region  is  called  a  mid-line 
segment  and  is  either  free  or  occupied.  A  mid-line  segment 
can  either  be  occupied  by  a  single  logic  node,  a  single  logic 
node  and  any  number  of  feedthrough  nodes,  or  any  number 
of  feedthrough  nodes.  The  region  is  bounded  on  the  left  and 
the  right  by  vertical  lines  called  net  separator  segments.  We 

^For  each  net,  netspan  —  pmax  —pminj  where  the  two  numbers  denote  the 
maximum  and  the  minimum  pin  level  of  the  net 
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(d)  Total  wire  crossing  (before  optimization)  (e)  Total  wire  crossing  (before  optimization)  (f)  Total  wire  crossing  (before  optimization) 


Fig.  2.  Evaluation  of  random  placements  of  V65E3-64. 


index  net  nodes  of  the  region  either  with  net  separator  seg¬ 
ment  i  —  1  or  net  separator  segment  each  being  of  height 

1.  The  region  jR(2,i)  is  bounded  on  the  top  and  the  bottom 
by  cell  separator  segments^  J  —  1  and  j.  We  call  the  left  half 
of  a  cell  separator  a  net-to-cell  separator  segment^  the  right 
half  a  cell-to-net  separator  segment  Each  has  a  width  of  0.5. 
A  cell  separator^  addressed  with  j  G  [OJmax],  is  a  union  of 
non-overlapping  cell  separator  segments.  A  net  separator^  ad¬ 
dressed  with  i  6  [Ojimaa:],  is  a  union  of  non-overlapping  net 
separator  segments.  We  may  also  refer  to  the  index  i  as  the 
level  of  net  nodes.  A  mid-line  is  a  union  of  non-overlapping 
mid-line  segments.  Each  mid-line  segment  is  addressed  by 
the  index  of  the  adjacent  net  separator  on  its  right  and  the 
adjacent  cell  separator  below. 

Elements  of  grid  region  i?(i,  j). 

1.  Two  cell  separator  segments,  shown  in  Figure  3  by 
dotted  lines.  One  is  at  the  boundary  of  grid  regions 
R{i,j  —  1)  and  R{i,3)  and  the  other  at  the  boundary  of 
grid  regions  R{ijj)  and  R{ij3  +  1). 

2.  An  array  of  net  nodes  of  type  B.  The  element  of  this 

array  is  labeled  (i  ”  1,  A;)b- 

3.  An  array  of  cell  nodes  and  feedthrough  nodes  of  type  A. 


The  element  of  this  array  is  labeled  {ijjjk)A- 

4.  An  array  of  net  nodes  of  type  B,  the  k^^  element  of 
which  is  labeled  (i, 

5.  Edges  connecting  the  nodes  in  the  region.  The  edges 

are  classified  as  type  A,  if  they  connect  nodes  labeled 
iiJiyki)A  and  {ij2,k2)B-  They  are  classified  as  type 
B,if  they  connect  nodes  labeled  {i  —  and 

{ijj2,k2)A- 

Properties  of  grid  region 

PI  In  each  grid  region  R{i,j)i  there  is  at  most  one  cell 
node. 

P2  For  a  cell  in  there  could  be  several  permissible 

orientations.  Each  orientation  defines  the  order  of  the 
pins  of  the  cell. 

P3  It  is  possible  to  derive  consistent  co-ordinates  for  each 
pin  of  a  node  using  its  label  and  orientation.  For  an  edge 
e  =  (“u,  t;),  we  define  y{u)  and  y{v)  as  the  co-ordinates  of 
the  pins  it  connects. 

P4  The  cell  separators  in  a  channel  induce  pseudo-edges. 
We  can  use  this  property  to  calculate  wiring  density. 
These  pseudo-edges  connect  nodes  which  are  at  the  in¬ 
tersections  of  net  separators  and  midlines.  The  co- 
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Fig.  3.  Elements  of  grid  region  j). 


The  I  factor  is  due  to  the  fact  that  each  crossing  is  counted 
exactly  twice,  once  for  each  of  the  two  edges  that  produce  it. 

Critical  wire  crossings  is  the  maximum  number  of  crossings 
that  appear  on  any  path  from  an  input  to  an  output. 

N-l 

wc.cri  =  wc(ui,ui^i)}  (7) 

i=l 


Wiring  Density  for  a  channel  C,  and  a  cell  separator  s,  is  the 
number  of  edges  that  cross  5  in  channel  C. 

wd{C,s)=^  C7’055(s,e)  (8) 

e^E(C) 

where  E{C)  is  the  set  of  edges  in  channel  C. 

Critical  wire  density  for  a  channel  C  (cell-to-net  or  net-to-cell) 
is  the  maximum  number  of  edges  that  cross  any  cell  separator 
in  that  channel. 


ordinates  for  these  nodes  are  derived  consistently  with 
those  of  other  nodes. 

P5  The  length  of  each  edge  e,  which  connects  nodes 
and  (i2,i2,fc2)B  is  defined  and  evaluated  as: 

j  /  \  /  0.5  if  ji  =  j2 

«en(e)-|  otherwise 

P6  Let  ei  =  (wi,vi)  and  e2  =  (u2,V2)  be  edges  such  that 
uifVi  are  in  the  same  wire  channel  and  are  pins  of  nodes 
with  the  same  subscript  {A  or  B)  as  U2  and  V2  respec¬ 
tively.  We  define  a  function  cross(ei,e2)  as  follows: 


wdjcri{C)  =  Max3es{'wd{CjS))  (9) 

Total  Wire  Density  for  a  placement  P  is  the  sum  of  the  critical 
wire  densities  for  all  the  channels. 

wdJtot{P)  =  wdjcri{C)  (10) 

C6C 

where  C  is  the  set  of  all  the  channels  in  the  placement. 

Average  Wire  Density  for  a  placement  P,  is  the  mean  of  the 
critical  wire  densities  for  all  the  channels. 


cross(ei,e2)  =  I  J 


if  (j/(tii)  -  y{tt2)){y{vi)  -  y{v2))  <  0 
otherwise 


(2) 

The  wire  crossing  number  for  each  edge  e,  is  defined  as 
follows: 


wc{e)  = 


E 


cro5s(e,ei) 


(3) 


eieT(e) 


awd{P)  = 


wdJot{P) 

|C1 


Width  of  a  channel  C  is  defined  as 

width{C)  =  1  +  0.5(tyd-cri(C)) 


(11) 


(12) 


where  T(e)  is  the  set  of  all  edges  of  the  same  type  as  e 
in  the  same  region  as  e. 

Figures  of  merit.  Let  E  be  the  set  of  edges. 

Total  wire  length  is  the  sum  of  the  wire  length  contribution 
of  each  edge. 

wlJtot  =  ien(e)  (4) 

Critical  wire  length  is  the  length  of  the  longest  path  starting 
at  any  primary  input  and  ending  at  a  primary  output.  This 
includes  feedback  primary  inputs  and  outputs  based  on  edges 
that  break  the  cycle. 

iV-l 

wljcri  =  Max(uiei,UNeo){^^  len{ui,  Ui+i)}  (5) 

*=i 


Width  of  the  placement  P  is  the  sum  of  the  widths  of  all  the 
channels. 

p.width{P)  =  width{C)  (13) 

cec 

Height  of  a  midline  M  is  defined  as 

height{M)  =  c(M)  +  0.1/(M)  (14) 

where  c(M)  is  the  number  of  cell  nodes  in  the  midline,  and 
f{M)  is  the  number  of  feedthrough  nodes  in  the  midline.  The 
0.1  factor  is  an  empirical  one. 

Height  of  a  placement  P  is  the  maximum  height  of  any  mid¬ 
line. 

pJieight{P)  =  MaxMeM{height{M)}  (15) 


where  I  is  the  set  of  primary  inputs,  O  is  the  set  of  primary 
outputs  and  ('Ui,Ui-f-i)  G  E. 

Total  wire  crossings  is  the  total  number  of  crossings  in  the 
placement. 

wcJtot  =  -  ^2,  ^c(e)  (6) 

eeE 


where  M  is  the  set  of  all  midlines  in  the  placement. 

j4 rea  of  a  placement  P  is  the  product  of  its  height  and  width. 

Area{P)  =  pJieight{P)  x  pjwidth{P)  (16) 

A  Placed  Netlist  Format.  To  describe  the  placement  of  a 
netlist  in  very  simple  terms,  we  found  it  convenient  to  extend, 
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.model  FA2 
.inputs  a  b  c 
.outputs  r  s 

#  wire  cut  0  b 

#  cell  cut  1  f 

#  wire  cut  1  d 

#  cell  cut  2  r 

#  wire  cut  2  r 

#  cell  cut  3  p 

#  wire  cut  3  s 
.names  h  d  r  # 
11  0 

#  wire  channel 

#  cell  channel 

.names  f  e  g  # 
11  0 

#  wire  channel 

#  cell  channel 


d  a 
d  e 

f  e  h  g  c 
g  h  c.lev2 
g  p  h  q  c.lev2 
s  q 

given  order 

hOl  dOl  r02 
hQ2  dOl  rQl 

reverse  order 

eQl  fQl  gQ2 
eQ2  fQl  g02 


.names  c  c.lev2  #  given  order 
1  1 

#  wire  channel  c01  c.lev2Q2 

#  cell  channel  c03  c,lev2Q3 

.end 


Fig.  4.  Illustrating  annotated  blif  format  for  Place-Eval. 

with  structured  comment  lines,  the  blif  format  [15].  The 
example  in  Figure  4  highlights  a  section  of  a  placed  blif  file 
that  corresponds  to  the  compacted  and  optimized  placement 
of  FA2,  full  adder  circuit. 

The  additional  comment  lines  #  wire  cut  i . . .  denote 
the  order  of  net  nodes  in  the  net  separator  i.  Similarly  the 
comment  lines  #  cell  cut  i . . .  denote  the  order  of  cell 
nodes  and  feedthrough  nodes  in  the  midline  i.  With  the  de¬ 
scription  of  any  cell  node  a  comment  of  #  given  order  or  # 
reverse  order  is  inserted  to  indicate  the  orientation  of  the 
cell.  Following  the  functional  description  of  each  node,  there 
are  2  comment  lines,  one  beginning  with  #  wire  channel 
which  gives  the  wire  channel  that  a  cell  node  and  its  inputs 
occupy,  the  other  begins  with  #  cell  channel  which  contains 
cell  channel  information  about  the  cell  node  and  its  inputs. 
Feedthrough  nodes  like  c .  Iev2  are  put  in  as  buffers. 

For  example,  the  line  #  wire  channel  eQl  fQl  gQ2  signi¬ 
fies  that  the  cell  node  g  is  located  in  the  wire  channel  2  and 
that  there  are  two  net  nodes  at  its  inputs,  e  and  f ,  which 
are  in  the  wire  channel  1.  The  line  #  cell  channel  eQ2  f  Q1 
g02  signifies  that  the  cell  node  g  is  located  in  the  cell  channel 
2  and  that  there  are  two  net  nodes  at  its  inputs,  e  and  f, 
which  are  in  the  cell  channels  2  and  1,  respectively. 

V.  TOCO  Placement  Algorithm 

The  acronym  TOCO  is  based  on  the  four  steps  of  the  pro¬ 
posed  placement  algorithm:  Topological  sort  that  places  net 
nodes  and  cell  nodes  on  the  unit-size  grid  in  level  order. 
Optimization  of  level  order  placement  for  niinimum  total  wire 
crossing  subject  to  the  level  order,  Compact  ion  and  folding 
of  the  level-order  placement  into  a  near-perfect  square  layout, 
Optimization  of  compacted  placement  for  minimum  total  wire 
crossing  subject  to  the  level  order  after  compaction,  and  min¬ 
imum  total  wire  density,  subject  to  order  for  minimum  total 
wire  crossing. 

To  minimize  the  total  wire  crossing,  we  currently  use  DOT, 
which  implements  a  variant  of  the  median  heuristic  [12].  To 
minimize  the  total  wire  density  after  total  wire  crossing  min¬ 


imization,  we  must  assign  net  nodes  and  feedthrough  cell 
nodes  to  grid  regions  such  that  critical  wire  density  is  min¬ 
imized  in  each  net-to-cell  and  cell-to-net  channel  -  subject 
to  node  order  imposed  by  the  total  wire  crossing  minimiza¬ 
tion  algorithm.  Our  current  approach  consists  of  making  a  k- 
way  partitioning  assignment  of  p  net  nodes  and  q  feedthrough 
nodes  to  k  grid  regions  in  each  channel,  such  that  the  wire 
crossing  on  the  kA  cell  channel  separator  segments  is  mini¬ 
mized,  subject  to  node  order  from  total  wire  crossing  mini¬ 
mization.  The  implementation  of  a  total  wire  density  opti¬ 
mization  algorithm  is  in  progress. 

Summary  Overview  of  TOCO. 

T  opological  sort.  This  procedure  may  include  the  secirch 
for  a  minimal  FVS  in  the  circuit  is  sequential.  It  is  the 
basis  for  generating  the  initial  placement  as  a  multi-level 
directed  bipartite  graph.  This  procedure  assigns  a  level 
to  every  cell  and  net  node  in  graph  G. 

O  ptimization  of  placement  P{G)  for  minimization  of  wire 
crossings.  At  this  point,  the  level  of  a  node  is  still  con¬ 
strained  by  its  level  in  the  topologically  sorted  graph. 

C  compact  and  fold.  Using  the  optimized  placement  ob¬ 
tained  from  the  wire-crossing  minimization  algorithm, 
we  fold  the  placement  into  a  compacted  placement.  Ba¬ 
sic  rules  are: 

R1  Remove  all  feedthroughs.  In  our  model, 
feedthroughs  have  no  function  other  than  that  of  rout¬ 
ing  nets  which  span  more  than  one  level. 

R2  Count  the  cell  nodes.  If  N  is  the  number  of  cell 
nodes,  then  the  number  of  wire  channels  in  the  com¬ 
pacted  form  is  The  number  of  cells  in  a  wire 

channel,  wire  channel  capacity ^  is  determined  simi¬ 
larly. 

R3  If  the  current  number  of  cell  nodes  in  wire  channel 
i  exceeds  its  wire  channel  capacity,  then  fold  forward 
(move  excess  cell  nodes  at  the  top  to  the  top  of  wire 
channel  i  -I- 1.) 

R4  If  the  current  number  of  cell  nodes  in  wire  channel 
i  is  less  than  its  wire  channel  capacity,  then  fold  back¬ 
wards  (move  deficient  cell  nodes  from  the  bottom  of 
the  next  available  wire  channel  to  the  bottom  of  wire 
channel  i.) 

R5  Insert  feedthroughs  %r  nets  that  span  more  than 
1-wire  channel. 

R6  Save  as  a  compacted  placement,  PJJompact{G). 

O  ptimization  of  placement  of  compacted  P{G).  This  con¬ 
sists  of  two  distinct  phases: 

PI  Minimize  the  total  wire  crossing  in  the  bipartite  form 
embedding  of  the  graph.  Net  nodes  and  cell  nodes 
are  constrained  to  the  levels  cissigned  to  them  by  the 
compaction  step. 

P2  Minimize  the  total  wire  density  in  the  bipartite  form 
embedding  of  the  graph.  Net  nodes  and  cell  nodes  are 
constrained  to  levels  and  the  order  assigned  to  them 
by  the  preceding  wire  crossing  minimization  phase. 
Wire  crossing  remains  unchanged  in  this  phase,  only 
the  total  wire  density  is  reduced. 

TOCO  implementation  and  Pseudocode.  Implementations 
of  various  phases  of  the  TOCO  procedure  have  been  optimized. 
Our  annotated  blif  format  seamlessly  introduces  net  nodes  for 
each  cell  node,  as  is  required  in  the  bipartite  form.  By  using 
the  same  names  for  cell  nodes  and  the  net  nodes  they  drive, 
the  complexity  of  any  phase  in  the  procedure  is  reduced  to 
that  of  searching  in  an  array.  Here  our  naming  convention 
allows  us  to  lexicographically  sort  arrays  to  make  multiple 
simultaneous  searches.  The  complexity  of  the  procedure  is 
dominated  by  these  sorting  calls  and  this  can  be  done  effi- 
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ciently  in  0{nlogn). 

Figure  5  shows  the  pseudo-code  for  the  TOCO  placement 
algorithm. 

procedure  TOCO-place(netlist) 

P(G)  =  Topological -Sort{nQt\\st) 

Pjopt{G)  =  Optimize{P(G)) 

PJCompact{G)  =  Compact{PX>pt{G)) 
Optimize{PX)ompact{G)) 

end 

procedure  Compact(netlist) 
Remove«Feedthroughs(netlist) 

N  =  number  JO  f  jcelljnodes 
number  JO  f  jwirejchannels  =  [>/N\ 
foreach  wire  channel 

if  cells  {wire -channel)  <  capacity  {wire  jchannel) 
FoldJBackward 
else 

FoldJForward 
Insert  JFeedthroughs 

end 

Fig.  5.  Pseudo  code  of  the  TOCO  placement  algorithm. 

VI.  Design  of  Experiments 

The  capability  to  synthesize  a  large  number  of  WSI  circuit 
mutants,  based  on  wire  perturbation  classes,  motivates  us  to 
examine  the  sampling  methods  that  arise  in  the  design  of  ex¬ 
periments.  Such  methods,  first  formalized  in  [16],  have  been 
adopted  widely  in  many  fields  of  science.  In  this  paper,  we 
adapt  them  to  analyze  the  performance  of  important  graph- 
based  algorithms  in  the  context  of  EDA.  For  each  reference 
circuit,  we  propose  to  synthesize  equivalence  subclasses  of  cir¬ 
cuit  mutants,  based  on  0  to  100%  perturbation.  Each  subclass 
contains  100  randomly  chosen  mutant  circuits,  each  listed  in 
a  different  random  order.  This  sample  size  is  large  enough  for 
the  sampling  distributions  to  be  considered  normal  or  nearly 
normal;  the  population  parameters  may  be  estimated  closely 
by  their  corresponding  sample  statistics.  The  eight  equiva¬ 
lence  subclasses,  labeled  from  A  to  jff,  are  defined  in  terms  of 
the  perturbations  we  use  to  generate  each  class.  In  order  to 
encourage  unbiased  experiments  with  these  classes,  we  have 
permuted  the  perturbations  relative  to  the  label  assignments: 

{A,B,C,D,E,F,G,H}  = 

permutation  (On;,  Iw,  2ty,  5%,  10%,  20%,  40%,  100%}  (17) 

In  (17),  WSI  classes  A-H  are  defined  either  in  terms  of  q- 
%  wire  perturbations  or  0-wire,  1-wire,  2-wire  {Ow,  Iw,  2w) 
perturbations.  Here,  we  only  disclose  that  the  0-wire  per¬ 
turbation  class  is  class _D.  We  plan  to  identify  the  la¬ 
bels  A’H  in  terms  of  the  respective  perturbation  classes  in 
(17)  later,  once  there  are  additional  experiments  reported 
by  others  and  participants  have  the  opportunity  to  meet 
and  present  their  results  at  a  joint  session  of  a  confer¬ 
ence.  More  details  about  such  plans  can  be  found  under 
http :  /  /  www.cbl .  ncsu .  edu  /  experiments  / . 

A  case  study  tutorial  of  average-case  performance  of  two  al¬ 
gorithms  and  their  differences  has  demonstrated  that  up  to  six 
distinct  equivalence  classes  of  data  are  useful  to  render  an  un¬ 
biased  comparison  of  two  well-known  sorting  algorithms  [17]. 
The  long-term  goal  of  this  series  of  experiments,  presently 
starting  with  eight  equivalence  classes,  is  to  facilitate  genera¬ 
tion  of  similar  comparisons  for  the  more  complex  and  diverse 
algorithms  in  EDA.  Whatever  may  be  decided  about  the  most 


(a)  algorithms-vs-mutants-vs-classes 

H-H_l  _  M-H_k  -  M-H>  I 


_  M-B_k 

_  M-B„b 

1 

M-A_l  . . 

. .  H-A_k  . . 

. .  M-A.b  1 

! 

1 

Alg_l-I 

A_ll-I 

. . . .  A_lk-I 

. . . .  A_lb-I  1 

1 

Alg_l~F 

A_ll-F 

_  A_lk-F 

_  A.lb-F  1 

1 

1 

1 

Alg_j-I 

A_jl-I 

. . . .  A_jk-I 

. . . .  A. jb-I  1 

1 

AlgJ-F 

A_jl-I 

. . . .  A.jk-F 

....  AJb-F  1 

1 

1 

1 

Alg_a-I 

A_al-I 

. . . .  A_ak-I 

. . . .  A^ab-I  1 

Alg.a-F 

A_al-F 

_  A^ak-F 

. . . .  A„ab-F  1 

(b)  classes-vs-mutants  (for  a  given  algorithm  Alg_j) 


Alg.j-F 

M-.^l  .... 

M-._k  .... 

M-._b 

M-A 

H-B 

AJl-F  ., 
B.jl-F  .. 

..  AJk-F  .. 
..  B_jk-F  .. 

. .  A, jb-F 
. .  B_ jb-F 

{H-A_j-F> 

<M-A^j-F> 

M-H 

H_jl-F  .. 

. .  jk-F  . . 

, . .  H_ jb-F 

<M-H^j-F> 

..<M-._jk-F>  ..{M-._jb-F} 


(c)  statistics  of  algor ithms-vs-mut ant  classes 


M-A 

M-B 

.  . .  M-H 

Alg_l-I 

•CM-A^1-I> 

•CH-B^1-I> 

..  <M-H_1-I} 

Alg.l-F 

{M-A_l-F} 

{M-B_l-F> 

..  <H-H„1-F} 

Alg^j-I 

{M-A,j-I> 

..  -CM-H.j-I} 

AlgJ-F 

■CM-A.j-F> 

<H-B_j-F> 

..  <M-H_j-F} 

Alg_a-I 

<M-A.a-I} 

<M-B^a-I} 

..  -CM-H^a-I} 

Alg_a-F 

-[M-A_a-F> 

<M-B^a-F} 

..  <M-H_a-F> 

Fig.  6.  Data  structures  and  classes  for  the  proposed  experiments. 

suitable  number  of  equivalence  classes  through  wider  partici¬ 
pation  later  on,  the  class  of  0-wire  perturbations  will  remain 
important.  As  demonstrated  in  this  paper  as  well  as  earlier 
[18],  the  objective  functions  used  in  a  number  of  graph-based 
algorithms  can  be  very  sensitive  to  the  order  of  nodes  in  the 
graph,  even  when  graphs  are  isomorphic.  In  our  experiments, 
the  100  net  lists  in  the  0-wire  perturbation  class  are  simply 
isomorphic  instances  of  the  reference  netlist  in  a  randomized 
order. 

Paraphrasing  the  context  of  the  traditional  treatments  and 
blocks  [19],  we  propose  to  archive  data  in  the  context  of  al¬ 
gorithms  and  equivalence  class  mutants  as  shown  in  Figure 
6(a).  For  each  of  the  ‘a’  algorithms  we  consider  ‘6’  mutants 
in  one  of  the  equivalence  classes  in  (17).  For  each  algorithm 
Alg_j  and  mutant  M-X_k,  X  €  {A, we  record  two 
observations:  the  initial  value  of  the  objective  function  tuple 
and  the  final  value  of  the  objective  function  tuple 
X_jk--F.  The  initial  value  corresponds  to  a  placebo  treatment 
of  the  mutant  M-X_k:  it  is  the  value  of  the  objective  function 
before  engaging  the  algorithm  to  optimize  it.  The  final  value 
corresponds  to  the  optimized  value  of  the  objective  function 
after  engaging  the  algorithm  to  optimize  it. 

A  number  of  analyses  can  be  performed  once  data  is 
archived  as  shown  in  Figure  6(a)  and  only  a  few  axe  dis¬ 
cussed  in  this  paper.  For  the  most  part,  we  shall  concentrate 
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on  analyzing  data  as  presented  in  Figure  6(b).  In  particu¬ 
lar,  for  samples  associated  with  each  algorithm  Alg_j  and 
mutant  class  M-X,  X  G  {A, . . . ,  if},  we  evaluate  the  95% 
confidence  interval  of  the  sample  mean,  the  sample  mean, 
and  the  sample  variances  as  tuples  {M“X_j-I}  and  {M“X_j-F} 
respectively.  We  summarize  such  evaluations  in  the  form 
shown  in  Figure  6(c).  The  next  section  provides  represen¬ 
tative  summaries  of  data  samples  we  generated  and  archived 
under  http://www.cbl.ncsu.edu/experiments/. 

VII.  Experiments 

This  section  summarizes  experiments  based  on  four  equiva¬ 
lence  classes  of  circuit  mutants  of  the  reference  circuit  C1355 
[20,  21].  The  four  classes,  with  100  mutant  circuits  in  each 
class,  are  based  on  wiring  signature-invariant  perturbations 
as  defined  in  Section  VI:  W-Class_B,  W-Class>C,  W_ClassJ), 
and  W-Class_E.  Notably,  W.Class  J)  is  the  0-wire  perturbation 
class  -  all  net  lists  in  this  class  are  isomorphic. 

Context  of  experiments.  The  context  in  which  we  propose 
to  perform  placement  experiments  is  shown  in  Figure  7.  The 
flow  may  be  executed  as  follows: 

1.  prepare  netlists  of  appropriate  equivalence  classes; 

2.  place  any  number  of  equivalence  class  netlists  with  any 
placer  e.g  TOCO; 

3.  archive  results  of  placements  in  a  common  grid-based 
format,  such  as  suggested  in  Figure  4; 

4.  process  grid-based  placements  by  either  or  both: 

(a)  a  common  placement  evaluator  (e.g.  place-eval  that 
can  process  results  of  any  grid-based  placement  and 
produce  a  standardized  report  in  terms  of  grid-based 
layout  parameters,  optionally  mapped  to  technology- 
specific  parameters  from  the  library  database; 

(b)  a  common  router  that  can  process  results  of  any  grid- 
based  placement  and  produce  a  layout  and  a  stan¬ 
dardized  report  in  terms  of  technology-specific  layout 
parameters. 

Ideally,  one  may  be  able  to  correlate  the  grid-based  param¬ 
eter  reported  by  place^eval,  such  as  total  and  critical  wire 
crossings,  total  and  critical  wire  length,  total  and  critical  wire 
density,  to  comparable  parameters  measured  in  actual  lay¬ 
outs.  If  the  correlations  are  reasonable,  the  grid-based  model 
may  be  calibrated  for  fast  estimation  of  cell-based  layouts. 

Finally,  execution  of  placement  experiments,  with  different 
algorithms,  as  suggested  in  Figure  7,  will  allow  us  to  study 
their  behavior  and  improve  them  further. 

Figure  8  —  Histograms.  Here  we  contrast  layout  parame¬ 
ters  reported  by  the  placement  and  routing  tool  OASIS  [13] 
and  grid-based  parameters  optimized  by  TOCO  and  reported 
by  place„eval  -  for  100  instances  of  the  mutant  circuits  in 
W_Class_B.  The  coefficient  of  variation  (mean/stdev)  for  lay¬ 
out  area  reported  by  OASIS  is  4.4%,  and  3.8%  for  the  min¬ 
imized  wire  crossing  reported  by  TOCO.  It  would  be  useful  to 
know  whether  the  placements  from  TOCO  would  improve  the 
layout  beyond  the  one  currently  shown  for  OASIS. 

It  is  clear  that  that  wirelength  and  total  layout  area,  as 
reported  by  OASIS,  are  closely  correlated.  Minimizing  the 
wirelength  does  minimize  the  area.  The  square  of  the  corre¬ 
lation  coeflScient  of  total  area  to  total  wirelength  as  reported 
by  OASIS  is  0.835.  However,  current  results  with  TOCO  show 
that  the  correlation  of  total  wire  length  and  total  wire  cross¬ 
ing  is  not  high.  We  attribute  this  behavior  to  the  sub-optimal 
assignment  of  net  nodes  to  the  respective  cell  channels.  We 
expect  the  correlation  to  improve  once  we  complete  the  im¬ 
plementation  of  the  channel  wire  density  minimization  algo¬ 
rithm. 

Figure  8  —  Table.  Here  we  summarize  results  of  experi- 


Fig.  7.  Proposed  context  of  placement  experiments. 


ments  with  4  classes  of  circuit  mutants  of  C1355.  Each  entry 
reports  the  95%  confidence  interval  for  the  mean,  the  sam¬ 
ple  mean,  and  the  sample  variance.  The  initial  values  refer 
to  values  of  parameters  after  compaction  with  TOCO.  The  fi¬ 
nal  values  are  after  the  wire  crossing  optimization  phase.  For 
all  these  classes,  the  procedure  dramatically  reduces  the  total 
wire  length  (TWL),  critical  wire  length  (CWL),  total  wire  cross¬ 
ing  (TWC)  and  critical  wire  crossing  (CWC).  However  the  im¬ 
provements  in  other  parameters  such  as  average  wire  density 
(AWD)  and  width  are  not  as  dramatic.  We  expect  to  report  im¬ 
proved  results  for  the  average  wire  density  (AWD)  and  width  in 
the  final  version  of  this  paper;  upon  completion  of  the  channel 
wire  density  minimization  algorithm  implementation. 

Statistical  interpretations.  The  reader  should  observe  the 
sensitivity  of  algorithms  when  evaluating  the  0-wire  perturba¬ 
tion  class  (W_Class_D)  in  the  table  of  Figure  8  and  elsewhere 
[5].  Experiments  implemented  in  this  paper  were  designed  to 
begin  addressing  issues  such  as 

•  Consider,  for  a  given  mutant  class,  (1)  sample  mean  and 
standard  deviation  of  the  unoptimized  objective  function, 
and  (2)  sample  mean  and  standard  deviation  of  the  ob¬ 
jective  function  optimized  via  algorithm  Alg_j.  We  have 
to  decide:  (HO)  is  the  difference  in  the  means  due  to 
chance,  or  (HI)  is  it  due  to  the  effect  of  algorithm  Alg^j? 

•  Consider,  for  a  given  algorithm  Alg^j,  (1)  sample  mean 
and  standard  deviation  of  the  optimized  objective  func¬ 
tion  in  terms  of  mutant  class  A,  and  (2)  sample  mean 
and  standard  deviation  of  the  optimized  objective  func¬ 
tion  in  terms  of  mutant  class  B,  We  have  to  decide:  (HO) 
is  the  difference  in  the  means  due  to  chance,  or  (HI)  is 
it  due  to  differences  of  the  two  mutant  classes? 

•  Consider,  for  a  given  mutant  class  (1)  sample  mean  and 
standard  deviation  of  the  objective  function  optimized  by 
algorithm  Alg_  j  1,  and  (2)  sample  mean  and  standard  de¬ 
viation  of  the  objective  function  optimized  by  algorithm 
Alg_j2.  We  have  to  decide:  (HO)  is  the  difference  in 
the  means  due  to  chance,  or  (HI)  is  it  due  to  different 
performances  of  the  algorithms  Alg^jl  and  Alg_j2? 

To  encourage  wider  participation,  complete  tables  of  all  data 
samples  summarized  in  this  paper  have  been  archived  on 
our  web  site  (http://www.cbl.ncsu.edu/experiments/).  The 
archives  are  updated  periodically  with  additional  experi¬ 
ments  and  cases  of  more  detailed  statistical  analyses  [14]. 
The  web  site  provides  an  open  forum  to  interested  re¬ 
searchers  for  further  sampling,  tests  of  significance  and  hy- 
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Experiments  with  Class  B  of  C1355 


OASIS  Wiralength 
(C) 


TOCO  Total  Wire  Crossings 


(d) 


W-ClassJB 

W-Class_C 

W-ClassJD 

W_Class-E 

TWL-I 

TWL-F 

(50004.8,  50668.2] 
50336.53,  1658.37 
[  9522.90,  9677.24] 
9600.07,  387.81 

[48548.4,  49173.2] 
48860.77,  1562.05 
[  9419.45,  9587.73] 
9503.59,  422.84 

[44596.5,  45163.9] 
44880.19,  1418.66 
[  8140.59,  8322.09] 
8231.34,  456.05 

[47207.6,  47814.1] 
47510.85,  1516.32 
[  9123.57,  9295.27] 
9209.42,  431.41 

CWL-I 

CWL-F 

[  1018.47,  1044.35] 
1031.41,  65.04 
[  266.49,  273.91] 
270.20,  18.67 

[  1024.21,  1046.71] 
1035.46,  56.52 
[  268.43,  275.53] 
271.98,  17.82 

[  1048.33,  1070.91] 
1059.62,  56.73 
[  246.02,  253.14] 
249.58,  17.87 

[  1038.85,  1062.91] 
1050.88,  60.46 
[  269.95,  277.23] 
273.59,  18.30 

TWC-I 

TWC-F 

[  37094.69,  37427.63]  [  36291.43,  36637.33]  [  33880.97,  34182.31]  [  35471.07,  35811.87] 

37261.16,  836.56  36464.38,  869,08  34031.64,  757.11  35641.47,  856.28 

(  9879.17,  10031.39]  [  9029.02,  9184.40]  [  2934.66,  3083.30]  [  7326.78,  7457.10] 

9955.28,  382.48  9106.71,  390.42  3008.98,  373.45  7391.94,  327.45 

CWC-I 

CWC-F 

[  1510.96,  1544.08] 
1527.52,  83.20 
[  611.70,  629.54] 
620.62,  44.82 

[  1538.95,  1569.53] 
1554.24,  76.82 
[  585.29,  604.89] 
595.09,  49.22 

[  1516.07,  1537.59] 
1526.83,  54.08 
[  279.20,  293.36] 
286.28,  35.59 

(  1553.45,  1582.39] 
1567.92,  72.71 
[  517.79,  533.43] 
525.61,  39.28 

AWD-I 

AWD-F 

[  13.28,  13,42] 
13.35,  0.36 
[  13.03,  13.29] 
13.16,  0.65 

[  13.30,  13.44] 
13.37,  0.33 
[  13.09,  13.33] 
13.21,  0.62 

[  13.84,  13.94] 
13.89,  0.27 
[  11.72,  12.02] 
11.87,  0.73 

[  13.50,  13.64] 
13.57,  0.34 
[  12.84,  13.14] 
12.99,  0.73 

Width-I 

Width-F 

[  314.09,  317.27] 
315.68,  7.97 
[  308.68,  314.33] 
311.50,  14.23 

[  314.61,  317.49] 
316.05,  7.22 
[  309.90,  315.30] 
312.60,  13.59 

[  326.36,  328.74] 
327.55,  5.98 
[  279.92,  286.30] 
283.11,  16.04 

[  318.99,  321.93] 
320.46,  7.39 
[  304.61,  310.99] 
307.80,  16.04 

Height-I 

Height-F 

[  35.62,  35.84] 
35.73,  0.56 
[  35.62,  35.84] 
35.73,  0.56 

[  35.64,  35.84] 
35.74,  0.51 
[  35.64,  35.84] 
35.74,  0.51 

[  35.27,  35.47] 
35.37,  0.50 
[  35.27,  35.47] 
35.37,  0.50 

[  35.67,  35.89] 
35.78,  0.57 
[  35.67,  35.89] 
35.78,  0.57 

Area-I 

Area-F 

[  11210.71,  11350.01]  [  11231.66,  11362.14]  [  11528.71,  11639.77]  [  11400.81,  11533.15] 

11280.36,  350.00  11296.90,  327.85  11584.24,  279.04  11466.98,  332.52 

[  11017.38,  11248.24]  [  11064.37,  11285.37]  [  9893.25,  10133.79]  [  10888.50,  11143.50] 

11132.81,  580.04  11174.87,  555.28  10013.52,  604.36  11016.00,  640.70 

Fig.  8.  Summaxy  of  placement  results  in  TOCO  for  mutant  classes  of  circuit  C1355. 


potheseSj  and  statistical  inference  of  existing  data  and  bench¬ 
marks,  as  well  as  for  contributing  new  cases  of  benchmarks, 


new  data  from  experiments,  and  new  cases  of  statistical 
analysis.  The  site  will  maintain  contributions  of  partici- 
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pants  either  as  hyperlinks  to  data  and  documents  on  par¬ 
ticipant’s  web  site,  or  new  archives  will  be  created  under 
http://www.cbl.ncsu.edu/experiments/. 

VIII.  Conclusions 

The  main  conjecture  in  this  paper,  that  good  column-based 
placements  can  be  produced  by  minimizing  two  wire  crossing 
numbers,  requires  access  to  an  open  routing  tool  that  will 
accept  a  user-defined  placement.  Ideally,  more  than  a  single 
routing  tool  should  be  used  in  such  a  study.  We  are  looking 
for  a  partner  to  assist  us  in  this  process.  The  large  number  of 
mutant  circuits  in  equivalence  classes  provides  a  sound  basis 
for  the  design  of  benchmarking  experiments  that  can  prove 
or  disprove  this  conjecture  at  a  sound  level  of  significance. 

Collaborative  web-based  experiments  are  being  initiated 
under  http://www.cbl.ncsu.edu/experiments/.  The  proposed 
placement  evaluator  place_eval  will  process  results  of  any 
grid-based  placement  and  produce  a  standardized  report  in 
terms  of  grid-based  layout  parameters.  We  hope  that  the  ap¬ 
proach,  after  testing  by  several  users  and  some  calibration, 
will  prove  useful. 

Acknowledgments.  We  appreciate  the  access  to  tools  used 
in  this  research:  DOT  from  AT&T  Research. 
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